Testing the Copernican Principle with Hubble Parameter 
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By way of expressing the Hubble expansion rate for the general Lemaitre-Tolman-Bondi (LTB) 
metric as a function of cosmic time, we test the scale on which the Copernican Principle holds in 
the context of a void model. By performing parameter estimation on the CGBH void model, we 
show the Hubble parameter data favors a void with characteristic radius of 2 ~ 3 Gpc. This brings 
the void model closer, but not yet enough, to harmony with observational indications given by the 
background kinetic Sunyaev-Zel'dovich effect and the normalization of near-infrared galaxy lumi- 
nosity function. However, the test of such void models may ultimately lie in the future detection 
of the discrepancy between longitudinal and transverse expansion rates, a touchstone of inhomoge- 
neous models. With the proliferation of observational Hubble parameter data and future large-scale 
structure observation, a definitive test could be performed on the question of cosmic homogeneity. 

PACS numbers: 98.65.Dx, 98.80.Es, 98.62.Py 



Introduction — The Copernican Principle (CP) is the 
hypothesis that we do not occupy a privileged position 
in the Universe. It leads to the Friedmann-Robertson- 
Walker (FRW) metric as the metric of the homogeneous 
and isotropic background spacetime However, one 
may not expect the CP to hold on all scales of cosmolog- 
ical interest, for both theory and observation shows that 
large-scale structure can emerge even if a homogeneous 
and isotropic initial background is assumed. Recently, 
the observed near-infrared luminosity function from a 
complete sample of galaxies indicates that the data can- 
not rule out the possibility of our vicinity being described 
by a void model I2| . In addition, the void model may also 
serve as a possible explanation to the emergence of ac- 
celerated expansion of the Universe without employing 
an exotic component dubbed 'dark energy'. To further 
investigate the ramifications of such a non-CP scenario 
and ascertain the possible existence of a local void, we 
shall consider a variety of other cosmological tests, as 
laid out in the following Letter. We mainly make use of 
the observational Hubble parameter data (OHD) which 
is independent of CMB and galaxy distribution measure- 
ment, and their observational properties have not been 
elucidated well enough in the inhomogeneous void model. 

LTB dynamics and the void model — The Lemaitre- 
Tolman-Bondi (LTB) line element reads 



ds^ = ~de 



A'{r,tf 
1 - k{r) 



dr^ + A^r,t)dn^ 



(1) 



where ' denotes d/dr, and fc(r) is associated with the 
spatial curvature. The Friedmann-Robertson- Walker 
(FRW) metric can be recovered by imposing A{r, t) — 
a(t)r and k{r) = kr^. 

From the LTB metric one can go on writing down and 
solving the dynamical equations for LTB void models. 



One notices along the way that the spherical symmetric 
configuration gives rise to two expansion rates 



A' 



(2) 



After choosing a gauge ^o('') = ''i ^-nd a homogeneous 
'bang time', one needs only the boundary conditions to 
finally obtain the evolution history(see Ref. 0, 0] for 
more detailed treatments). Expressed as two functions, 
rim and -ff^Oj these boundary conditions define an LTB 
void model. Throughout this work, we employ the Con- 
strained GBH (CGBH) model i, in which 



an(r) = l + (f^o-l) 



1 - tanh[(r - ro)/2Ar] 
1 + tanh(ro/2Ar) 



(3) 



where fig describes the density at the symmetric center. 
To is the characteristic size of the void, and Ar describes 
the steepness of the void near the edge. 

To illustrate how the universe and its evolution in 
CGBH model look like, we choose VIq = 0.05, Hq = 
73 km s^^ Mpc^^, ro = 6 Gpc, Ar = 0.1 ro, and plot in 
Fig. [T] and Fig. [2] the density profile on different cosmic 
time (in Gyr) slices and on the light cone, respectively, 
and in Fig. [3] and Fig.|4]the profiles of the two expansion 
parameters on time slices and the light cone, respectively. 

In a homogeneous universe one always has iJy = H^. 
Observation of a difference, like that shown in Fig. 21 at 
any redshift, would imply spatial inhomogeneity. 

Observational Hubble parameter data — There are 
three main methods to measure the Hubble parame- 
ter H{z): by measuring the differential age of passively 
evolving galaxies (differential age method) [6l4llj , by mea- 
suring the baryon acoustic oscillation (BAO) along the 
line-of-sight direction from the spectroscopic galaxy sam- 
ple 1^, and the dipole of the luminosity distance of 
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gravitational wave sources (luminosity dipole method of 
standard sirens) 13, The radial BAO size method 
depends on the detailed evolution of perturbations not 
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FIG. 1. Density profiles at difTerent cosmic times. Time is in 
unit of Gyrs, and tlie densities are normalized to the back- 
ground density. It is seen that the void gets deeper and its 
shell gets denser as the universe evolves. 




1+z 



FIG. 2. Densities on the past light cone, normalized to the 
value at the void center. 

well understood in the LTB cosmology, although pro- 
gresses have been made [H, The luminosity dipole 
method of standard sirens by now has produced no obser- 
vational data yet. Therefore, the OHD used in this work 
refer exclusively to that by the differential age method. 

The Hubble parameter for FRW models with scale fac- 
tor a reads 

H ^ ^ '^^ (4) 

a I + z dtct ' 

where dtct is the variation of the cosmic time due to a 
small change in the redshift dz. For any galaxy one has 
Tca{z) = Tp + TcAiz), which simply states that the 



FIG. 3. Comparison of the H± and profiles at difTerent 
cosmic times. Time is in unit of Gyr, and the expansion 
parameters are normalized to the value at the void center. 
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FIG. 4. Comparison of the H± and i^y at difTerent redshifts 
(on the past light cone) . 

cosmic age Tca at redshift z equates the summation of 
the formation time of the galaxy, Tp, and the age of this 
galaxy, Tqa that can be determined spectroscopically. If 
we could find a group of galaxies that share a uniform for- 
mation time, i. e. Tp = const., we would then get a han- 
dle of dtct by simply measuring the age difference of those 
galaxies: dtct{z) = dTcA = dTp + dTcAiz) = dTcA^z). 
The passively evolving galaxies can be identified by fig- 
uring out at every redshift the oldest galaxies, which to- 
gether define the 'red envelop'. One assumes in this pro- 
cess the oldest galaxies formed at the same time (stan- 
dard cosmic chronometers), which is a natural assump- 
tion in an FRW universe (one may call it the galaxy- 
formation version of the cosmic Copernican Principle). 
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Of the two expansion rates defined in Eq. ([2]), the lon- 
gitudinal expansion rate H\\ turns out to have the same 
form as Eq. @ H\\ = — [l/{l-\-z)\{dz/dtct), hence corre- 
sponds to the observed H{z) 

The problem of using OHD in LTB models is that the 
basic assumption that the oldest galaxies share a same 
formation time might not hold any more, as discussed re- 
cently in Ref. because the background in LTB mod- 
els has considerable inhomogeneities. However, we argue 
that OHD is still valid in our context (see Discussion). 
In the following we will use the latest 23 data entries as 
listed in Refs. 

Constraints on the void model — We perform the 
constraints on the void model by both OHD and 
the background inhomogeneity-induced kSZ (BIkSZ) 
effect 171, |l8|. The direction dependent, therefore observ- 
able, temperature shift reads 

ATbi = TcMB X r ^e{n,z) MllillZldre , (5) 



where Tcmb = 2.73 K, Zg = 100 (the result essentially 
does not change as long as z(ro) ^ Ze), and 

vh{z) « (r(z), t{z)) - {r{z,),t{z))]A{r{z),t{z)) 

(6) 

We calculate BIkSZ power spectrum AT|j at I = 3000 
and its constraints on the (ro, fio) parameter plane, with 
Ar/ro = 0.21, Hq = 74 km s"^ Mpc"^ fixed at their re- 
spective best-fit values, which are in turn obtained from 
the Hubble parameter dataset. The resulting contours, 
as well as confidence regions from the OHD and the su- 
pernovae Union2 dataset, are plotted in Fig. [SJ The 
OHD data favor a smaller (and more tightly constrained) 
void than what the supernovae Union2 data do. Indeed, 
there is a clear discrepancy between those two datasets, 
as found in Ref. Q . This is a sign of inadequacy for this 
specific LTB model. Also, as pointed out in Ref. [l^, one 
can see from Fig. [5] that the Gpc-sized voids, as those fa- 
vored by the supernovae data, are incompatible with the 
BIkSZ measurement, hence are largely excluded. Now 
we can tell from the figure that the OHD dataset give 
slightly weaker, though basically the same, conclusion. 
Furthermore, the observed normalization of the near-IR 
galaxy luminosity function indicates that a void, if ex- 
ists, amounts to a few hundred MpcsQ- This could in 
principle be consistent with the BIkSZ measurement. 

Future BAO constraint — As is shown above, unlike 
the homogeneous cosmological models, H\\ can differ from 
i/^in LTB models. Therefore one straightforward way is 
to define the ratio between these two expansion rates £ — 
H\\/H±, which always equals to 1 for the homogeneous 
cosmological models, but deviates from 1 in LTB models. 
Specifically, we can further write £ as 
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FIG. 5. Background inhomogeneity-induced kSZ effect ATgi 
in fiK^ . The line marked with 7 corresponds to the 95% 
observational upper limit of the South Pole Telesco pe [l9l . |20| ] 
which is significant to the constraint of void models [131 • Also 
plotted are the la and 2a confidence regions from the OHD 
as well as those from the supernovae Union2 dataset. The 
vertical (red) bar marks the characteristic void radius r — 
250 Mpc in the small-void models of Refs. [U, [23] that are 
comparable with recent observational data of near-infrared 
galaxy luminosity function Q]. 



A H'^ 

'a'IiZ' 



(7) 



since H\\ = ^^H^ + ^H'^. 

Therefore, the violation oi £ — 1 could also be an in- 
dicator of LTB-type models. To adopt this criteria how- 
ever, one needs independent measurements of and 
H± at the same redshift or just the variation of H± on 
the light cone. The BAO feature imprinted in the non- 
relativistic matter such as galaxies distribution yields a 
further geometric test of homogeneity, and future large- 
volume BAO surveys will also allow us detect the BAO 
scale in both radial and transverse directions. Thus, we 
expect future BAO measurement would supply the in- 
formation of the criteria £ by radial Hubble parameter 
and transverse Hubble parameter H±, and improve 
greatly the testing of violation of homogeneity. 

Discussion — In our calculation of OHD constraints 
on the void model, we assume that OHD could be used 
in LTB models. Actually, the same formation time of 
the oldest galaxies is a basic assumption in obtaining 
OHD. However, in LTB models where the universe has 
a considerable background inhomogeneity, this assump- 
tion becomes unreasonable and some arguments are also 
recently given in Ref. [2^ where galaxy ages are used. 
First, despite the overall uniform- formation- age assump- 
tion, the validity of an H(z) data point requires the same 
formation time only inside the redshift bin where OHD is 
locally defined and obtained (Eq. [J] ) , even if the global 
density - hence the formation time of the oldest galaxies 
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- at different redshifts varies much. Secondly, a stan- 
dard viewpoint (referred to as the onion approximation) 
is to treat the LTB void universe as a group of thin shells 
structured together, and inside each of these spherical 
shells the matter is homogeneously distributed [2J|. For 
the OHD used in this paper, the size of each redshift bin 
is between 0.1 and 0.15, where the first limit is so cho- 
sen that the age evolution between the two bins is larger 
than the error in the age determination As the pre- 
cision of the age determination improves, we expect an 
even smaller bin size. To be sure about the validity of 
OHD used in LTB models, one needs the exact knowl- 
edge about the thickness of the shell given the size of a 
redshift bin, as well as the steepness of the density profile 
at the time the oldest galaxies formed. We discussed this 
issue in Ref. 4*1 . 

On the other hand, future observation is expected to 
yield ^ 2000 measurements for passively evolving galax- 
ies in the redshift range < z < 1.5 in the future 0. 
It has been estimated that about 1000 OHD entries at 
a 15% accuracy level will be determined with 10% error 
of the galaxies ages. In Ref. the power of OHD in 
the context of ACDM model has been assessed. With 
the increase of high quantity OHD, the power of OHD 
constraining void model should also be greatly improved 
for constraining the void models 

Although the large void model appear to be ruled out 
by some cosmological observations, future OHD measure- 
ment in both radial and transverse directions, as an alter- 
native and complementary cosmological test, could give 
a tight constraint on LTB model with a small void. If the 
transverse BAO information can be realized from future 
large-scale structure observations, we should be able to 
arrive at a definite test of spatial homogeneity of the Uni- 
verse. In this context, the role played by the transverse 
BAO is complementary to the radial BAO discussed in 
Ref. fi^. 
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